Staffing many- server systems 
with admission control and retrials 



A.J.E.M. Janssen* J.S.H. van Leeuwaarden^ 

February 14, 2013 

Abstract 

In many-server systems it is crucial to staff the right number of servers so that targeted 
service levels are met. These staffing problems typically lead to constraint satisfaction problems 
that are hard to solve. During the last decade, a powerful many-server asymptotic theory 
has been developed to solve such problems and optimal staffing rules are known to obey the 
square-root staffing principle. This paper develops many-server asymptotics in the so-called 
QED regime, and presents refinements to many-server asymptotics and square-root staffing for 
a Markovian queucing model with admission control and retrials. 
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1 Introduction 

A key challenge in large many-server systems is to balance the trade-off between operational costs 
and quality-of-service offered to customers. There is by now a vast literature on the asymptotic 
analysis of many-server systems. Many papers describe limiting behavior of a sequence of queues, 
which is then used to approximate the characteristics of a member of the sequence, i.e., the per- 
formance of a finite-sized queueing system. Depending on how this sequence is parameterized, its 
limiting behavior is different, giving rise to different approximations (see [3] and [12]). One of 
the most popular approximations arises in the Quality-and-Efficiency-Driven (QED) regime, also 
known as the Halfin-Whitt regime [8] , in which the number of servers s and the offered workload A 
are related according to a square-root principle s = A + /3V% for a constant /3, and s (and A) are 
taken to infinity. 

In this paper we consider a Markovian many-server system with admission control in the QED 
regime. We consider an admission control policy that lets an arriving customer enter the system 
according to a probability depending on the queue length. In particular, a customer meeting upon 
arrival k other customers is admitted with probability pk, and we shall allow for a wide range 
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of such admission policies characterized by (pk)k^s- For such admission-controlled systems, we 
consider the QED regime and the square-root staffing rule A = s — jy/s for some fixed constant 
7. This is asymptotically (as s — > 00), and for all practical purposes, equivalent with setting 
s = A + PVX for some fixed constant (3. The reason for staffing via A instead of s is that it allows 
for mathematically more elegant derivations. All results in this paper regarding A = s — "fy/s can 
be converted into results for s = A + (3y/\ at the cost of increased notational burden. 

Since certain customers are rejected upon arrival, it seems natural to extend the model with 
the feature that rejected customers can reattempt. The modeling of reattempts or retrials is known 
to be challenging [4, 6], which is why one often resorts to computational approaches [1]. These 
numerical approaches face increasing numerical difficulties when the number of servers becomes 
large, which is precisely the regime we are interested in. Therefore, we combine the QED regime 
with a limiting regime for slow retrials, meaning that rejected customers reattempt after a relatively 
(compared to the time scale of the system) long time. The combination of these two asymptotic 
regimes leads to a tractable model with a closed-form solution, for which we are able to derive QED 
approximations for some of the relevant performance measures. 

We leverage these QED approximations to obtain results for staffing problems. The QED 
regime is particulary suited for staffing many-server systems, as it combines large capacity with 
high utilization. This appealing feature was exploited in [3] that introduced a rigorous asymptotic 
framework for applying the square-root staffing principle to two classes of problems: constraint 
satisfaction and cost minimization. In [3] it was observed that square-root staffing based on QED 
approximations is effective and accurate over a wide range of system parameters for the Erlang 
C model. We apply the same approach to develop square-root staffing rules for several constraint 
satisfaction problems for the system at hand. 

Square-root staffing. The core of staffing problems in many-server systems is to determine the 
right trade-off between quality and capacity. Quality is formulated in terms of some targeted service 
level. Take as an example the delay probability Dp(s,X) (see (2.8)). A large delay probability is 
perceived as negative, and the targeted service level could be to keep the delay probability below 
some value s. The smaller e, the better the offered service. Once the targeted service level is set, 
the objective from the system's perspective is to determine the largest load A (or the lowest staffing 
level s) such that the target Dp(s, A) ^ e is met. This is what we have referred to as a constraint 
satisfaction problem. 

The delay probability is a function of the two model parameters s and A, and of the admission 
policy. Because the delay probability is a continuous and monotone increasing function in A, the 
constraint satisfaction problem is equivalent to finding the A op t such that Dp(s, A op t) = £• To solve 
this inverse problem, we shall invoke the theory of asymptotic dimensioning introduced in [3] and 
extended in the first part of this paper to admission control and retrials. This theory fully exploits 
the QED regime for large systems, in a way that reduces considerably the complexity of the inverse 
problem. That is, in the QED regime (with A = s — 'yy/s, s — > 00, 7 not scaling with s), the perfor- 
mance measures in our model can be approximated by their limiting counterparts. For instance, 
Df(s,\) can be approximated by some function -D*(7) that only depends on 7 (and no longer on 
s or A). Hence, the inverse problem can then be approximatively solved by searching for the 7* 
such that D*(7*) = e, and then setting the load according to A* = s — j^y/s. We refer to this rule 
as conventional square-root staffing. In this asymptotic approach, one expects that the better the 
approximation Dp(s, A) « ^(7), the smaller the error |A opt — A*|. Based on the QED regime, one 
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also expects the approximation A* to be accurate for large values of s, i.e., for large-scale systems 
with many servers. Indeed, we prove that |A opt — A*| = 0(1), where a function /(A) = 0(g(X)) if 
limsup A _ ) . 00 \f(X)/g(X)\ < oo. 

Refined staffing. We shall also derive refined staffing rules, for which we first develop refined 
QED approximations for the objective function, and then characterize the approximative solutions 
to the constraint satisfaction problems. The refined staffing rules are of the form 

A. = s - ^yfs + r., (1.1) 

with r, a simple function of 7* , s and e. We shall uniquely identify r, (for each considered constraint 
satisfaction problem), and prove that the refined staffing level in (1.1) satisfies 

Aopt - A. = 0(s-^ 2 ), (1.2) 

We refer to the order term that expresses the difference between the exact optimal staffing level and 
the approximate staffing level as the optimality gap. Hence, the optimality gap of A, is 0(s~ 1 ^ 2 ), 
which suggests that the staffing level A, becomes more accurate as s increases. Note that A. = 
A* + 7». Since the optimality gap of the conventional staffing level A* equals O(l), we can expect 
that A, should be a more accurate prescription than A*. In addition, because 7, in fact describes 
the optimality gap of A* , or more precisely, 

Aopt -A, =7. + 0(^r 1/2 ), (1.3) 

it allows us to perform an analytical assessment of the accuracy of both conventional and refined 
square-root staffing. 

The influence of admission policies and retrials. Our main results are captured in Theorems 
7 and 8, which formally establish the necessary results for both conventional staffing rules and 
their refinements. These two theorems establish that as the system size s increases, the difference 
between the true maximal load that adheres to the constraint A op t, and the conventional square- 
root staffing prescription A*, converges to the real number 7., see (1.3). We stress that these results 
hold across the large class of admission control policies considered in this paper, and for the systems 
with and without retrials. 

We obtain the explicit expression of 7, in each case, and show that it is always a simple function 
of 7* . We can therefore investigate numerically the impact of the refinements in comparison with 
conventional square-root staffing. It turns out that, unlike in the classical Erlang C model [3], the 
refinement 7, is significant in many cases, due to different system parameters. Our findings further 
suggest that in the presence of admission control, more care needs to be taken in applying many- 
server asymptotic results to small or moderate size systems. In particular, when the admission 
control becomes more lenient, and hence on average more customers are admitted, our numerical 
findings suggest that the refinements to the conventional square-root staffing are necessary. 

Overview of the paper. In this paper we apply the QED regime to a Markovian many-server 
system with admission control and retrials. We analyze this system in the limiting regime, in a 
similar spirit as was done for the Erlang C model [8, 11], the Erlang B model [10] and the Erlang A 
model (with abandonments) [14]. In the recent paper [2] on a many-server system with retrials, the 
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admission control policy that rejects all delayed customers (loss system) was analyzed. The current 
paper strongly builds on several results obtained in [2]. Compared to earlier studies [8, 11, 10, 14, 2], 
the system in this paper brings about additional mathematical challenges, because of the effects of 
rejection and reattempts. In short, we make the following contributions: 

(i) We consider two stationary performance measures: The probability that an arriving cus- 
tomer finds all servers occupied Dp(s,X), and the probability that an arriving customers is 
rejected Dp(s, A). For both performance measures we derive in Theorems 3 and 4 the limiting 
expressions in the QED regime, for the cases with and without retrials. 

(ii) We next consider dimensioning problems of the type: For fixed s, find A for which the rejection 
probability Dp(s, A) has some prescribed value. We solve this inverse problem both for finite s 
and for s — > oo in the QED regime. The inverse problem in the QED regime is easier to solve, 
and provides a leading-order estimate, in terms of the hazard rate of the standard normal 
distribution, for the finite-s inverse problem. We characterize explicitly the error made by 
replacing the finite- s inverse problem by its QED counterpart, by deriving correction terms in 
Theorems 7 and 8. These correction terms are then also used to develop the refined staffing 
rules and establish the optimality gaps. 

(iii) We then consider dimensioning problems of the type: For fixed s, find A for which the 
carried traffic A(l — Dp(s, A)) has a prescribed value. These dimensioning problem are shown 
to have typically two solutions: A small A leads to a large probability of being admitted 
1 — Dp(s,X), while a large A leads to a small admittance probability; see Theorems 9 and 
11. This phenomenon of two solutions has also been observed in certain loss networks with 
alternative two-link routing [7]. 

(iv) The more technical part of this paper deals with incorporating the effect of retrials, for which 
an essential role is played by the generalized Cohen equation Q = (A + Q)Dp(s, A + 0) and 
its solution representing the average rate of retrials. In the next section we explain in more 
detail the fundamental role of this equation in our study. 

The paper is structured as follows. In Section 2 we describe the model. In Section 3 we present 
the limiting expressions for the performance measures in the QED regime, and we use these results 
to deal with the dimensioning problems in Section 4. Finally, we present in Section 5 results for 
dimensioning problems with multiple solutions. The proofs of the main results are given in Sections 
6-8, while the proofs of supporting results are deferred to the appendix. 

2 Model description 

The basic model described in this section is taken from [5]. Consider a system with s parallel 
servers to which customers arrive according to a Poisson process with rate A. The service times of 
customers are exponentially distributed with mean one. An admission policy dictates whether a 
customer is admitted to the system or rejected. A customer that finds upon arrival a free server is 
immediately assigned to that server, and leaves the system after service. A customer that finds upon 
arrival k other customers in the system, k ^ s, is allowed to join the queue with probability p^ and 
is rejected with probability 1 — Pk- In this way, the sequence (pk)k^s defines the admission policy. 
Since we are interested in large many-server systems working at critical load, and hence serving 



4 



many customers, the parameter pk should be interpreted as the fraction of customers admitted, 
instead of the probability that determines the fate of just one single customer. For the results 
presented in this paper we impose only mild conditions on the sequence (pk)k^s, allowing for a 
wide range of admission policies to be considered. 

Under the above Markovian assumptions, and assuming that all inter arrival times and service 
times are mutually independent, the system can be described as a birth-and-death process (C(t))t^o 
with C(t) the number of customers in the system at time t. The birth and death rates from state 
k are Xpk (with pk = 1 for k = 0, . . . , s — 1) and min(A:, s), respectively. Assuming the stationary 
distribution exists, with ir^ = lim^oo P (C(i) = k), it readily follows from solving the detailed 
balance equations that 



with 



k=0 



and 



F{x) = Y,Ps---Ps +n x n+1 . (2.3) 

n=0 

From (2.2) it can be seen that the stationary distribution exists when F(^) < oo. Since p s + n G [0, 1], 
we have that F(j) < oo when A < s. When A ^ s we need to be more careful. The radius of 
convergence of the power series F(x) is given by 1/P with 
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P := lim sup (p s ■ ... • p s+n ) n+i g [0,1]. (2.4) 



In a major portion of the main text, we assume the following condition: 

P G [0, 1), F(— — 0) = lim F(x) = oo. (2.5) 
p ast^ 

Under this condition it can be easily observed from (2.3) that the stability condition for our system 
becomes 

A G [0,A P ) with A P = p (2.6) 

where Ap = oo when P = 0. The condition (2.5) is certainly not as general as possible to develop the 
theory, but it excludes cases that need separate consideration, thereby distracting attention from 
the bottom line of the exposition. Stability is guaranteed when A < s. Also, when limk^ooPk = 0) 
we have P = and thus stability for all A 0. Condition (2.5) is also satisfied for the case 
Pk = p G (0, 1) for all k ^ s, where ^(x) = px/(l —px) so that -^(-p — 0) = oo, with P = p, indeed. 
We exclude at this point the case P = 1, which would for example occur in the case pk = 1, ^ s, 
and the cases that F(p — 0) < oo. However, in Section 6 and Appendix A, the results are proved 
under general conditions. 
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Let 

denote the Erlang B formula representing the stationary blocking probability in an M/M/s/s 
system. A crucial performance measure is the stationary probability Df(s,X) = X^a^= s n k that an 
arriving customer finds all servers occupied, given by 

l + F( 



D-As,» = - (2.8) 



where Dp 1 is short-hand notation for (Dp) -1 . Note that for p^ = 0, k ^ s, the term F(X) vanishes, 
and the probability Dp(s,X) reduces to B(s,X). Also, for p^ = 1, k ^ s, the probability Dp(s,X) 
reduces to the Erlang C formula given by 

C(s,X) = 1 pfa , (2.9) 

2^k=0 A / Kl + ( s _i)!( s _A) 

representing the stationary delay probability in an M/M/s system. 

Another important performance measure is the stationary probability Dp(s,X) = ^2'kL s 'KkPk 
of being rejected, given by 

B-Hs,X) + F(i 
1 + (1 — s/X)F(j) 

It follows from results in Section 6 that 



max{0, 1 - -} ^ D§ (s, X) ^ B(s, X) ^ D F (s, X) ^ 1 (2.11) 
A 

for ^ A < Xp. Moreover, as a consequence of condition (2.5) we have 

D F (s,X P -0) = 1, Dp(s, Xp — 0) = 1 — P. (2.12) 

The birth-and-death process (C(t))t^o relies on the assumption that rejected customers are 
considered lost. Alternatively, we could assume that rejected customers reattempt to enter the 
system after some time. In that case, rejected customers start producing reattempts until they are 
allowed to enter. Assume that periods between successive reattempts of a rejected customer are 
exponentially distributed with rate fi, independent of interarrival and service times. The system can 
then be described as a two-dimensional process (C(t), N(t))t^o with C(t) the number of customers 
in the system and N(t) the number of rejected customers at time t. Under the above assumptions 
this process is a continuous-time Markov chain on the lattice infinite strip {0, 1, . . . , s} x Z+. 

Since the transition rates of this process clearly depend on the second coordinate, the process 
{(C(t), N(t));t ^ 0} is hard to analyze. In fact, even deriving the stationary distribution poses 
analytical difficulties, and no closed- form solution seems to be available. We therefore make the 
following assumption: Reattempts arrive to the system according to a Poisson process with rate 
$7, independent of the Poisson process of customers that arrive to the system for the first time. 
This assumption is also known as the retrials see time averages (RTA) approximation. Under this 
assumption, the total flow of customers arriving to the system is a Poisson process with rate A + 0. 
The unknown rate f2 should then be the solution to 

n = (x + n)D^( s ,x + n). (2.13) 
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Equation (2.13) is intuitively clear as it equates two expressions for the rate of reattempts. In some 
cases the RTA approximation can be theoretically justified. Cohen [4] showed that the system 
with pk = for k s, in the limit as fx 4- 0, behaves as an Erlang loss system, except with an 
increased arrival intensity. More specifically, as [i J, 0, the distribution of the number of busy servers 
converges to the corresponding distribution for the standard Erlang loss system M/M/s/s (which 
is a truncated Poisson distribution), but with increased arrival rate A + J7, where is defined as 
the solution to (2.13). Indeed, in the case of infinitely long retrial times, it fits intuition that the 
flow of reattempts is independent from the flow of primary customers. For retrial queues with finite 
retrial times, the RTA approximation has proved useful and accurate for many retrial systems. We 
shall refer to (2.13) as the generalized Cohen equation. 

Theorem 1 (Unique solution Cohen's equation). Under condition (2.5), there is a unique solution 
^s,f(i) °f equation (2.13) for any A 6 (0, s). 

A proof of Theorem 1 can be distilled from [5, Section 3]. In Appendix A we shall present a 
self-contained proof of Theorem 1 under conditions more general than (2.13). 



3 QED limits 

The QED regime for many-server systems refers to scaling of the arrival rate A and the number 
of servers s such that, while both A and s increase toward infinity, the traffic intensity p = X/s 
approaches unity and 

(l-p)v^ 7 , (3.1) 

where 7 is a fixed constant. The scaling combines large capacity with high utilization. For the 
Erlang loss and delay systems, this kind of scaling leads to the classical results (see e.g. [10, Section 
5.2]), for 7 G (— oo,y/s) fixed, 

1 1 $(7) / 1 \ 

— B- 1 (s,s- 1 V~s) = ^r J : + (-= , s^oo, (3.2) 
V s y(7J v Vs/ 

and for 7 £ (0, y/s) fixed, 

lnnC( S , S -7v^)=(l+7|^)~ 1 , (3.3) 

where $>(x) and y>{x) denote the standard normal cumulative distribution function and density, 
respectively. 

The following result will prove useful in establishing QED limiting results. 
Lemma 2 (Two decompositions). For A G [0, Xp), 

D-\s, A) = (1 - q x )B- 1 (s, A) + q x C-\s } A), (3.4) 

D- R (s, A) = B-\s, A) + C- 1 ^, A), (3.5) 

1 - Qx 

where 

Qx = X s { ■ (3-6) 
1 + ^(1) 
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The proof of Lemma 2 is given in Appendix A. It can also be obtained from more general results 
in [13]. Note that the function C(s, A) is also defined for A > s, while in the M/M/s queue the 
stability condition is A < s. It is further shown in Appendix A that q\ ^ 1, with q\ = 1 if and 
only if pk = 1, k ^ s. Thus, for instance, Dp R always exceeds -B -1 and the excess is given by the 
second term of the right-hand side of (3.5), which is the product of a factor entirely determined 
by the admission policy and the Erlang C formula. Also, Dp 1 is a convex combination, with a 
7-dependent convexity parameter 1 — q\, of the Erlang B and C formulae. When an admission 
policy is mild, implying that q\ is close to 1, we have that the Erlang C formula is dominant. 
When an admission policy is strict, the Erlang B formula is dominant. Aside from these general 
comments, the variety of weight functions q\ that can occur, see (3.6), is rather substantial. In the 
QED regime, though, the Erlang B formula is always dominant. Indeed, in the QED limit, we have 
X/s —7- 1, and so q\ —> F(l)/(1 + F(l)) which is a finite number G (0, 1) by our condition (2.5). 
Now (3.2) and (3.3) show that B~ l grows like y/s while C~ l remains bounded as s — > oo. 

We now apply the scaling (3.1) to the system with admission policy. We henceforth keep working 
with the notation for the QED regime in (3.1), which is why we reformulate the stability condition 
(2.6) as 

7G(7P,v^] (3.7) 

with 

7P = -^=^^e (-oo,0). (3.8) 

Theorem 3 (QED limits without retrials). Under condition (2.5), for A = s — 7\/i, with 7 £ 
(—00, y/s\ fixed, 

hm v^M*,«-7V^ = (i + ^(i))lrT' ( 3 - 9 ) 

hm ^D§(s, s -^) = ^r\- (3-10) 
Proof. We have by continuity of F(x) at x = 1 that for fixed 7 E (—00, y/s\ 

f(^^)=f(i-^L)=F(1) + o(1), s^oo. (3.11) 

1 1 1 $(V) 

D F \ S ,s- 1 yfs) = — — - -All + (1), s^oo, (3.12) 
1 + F(1) y>(7) 

J-D F %,s- 1 ^) = ^ + o(-^), s^oo. (3.13) 
This implies (3.9) and (3.10). □ 

A first observation is that the limiting expressions (3.9) and (3.10) are similar as for the Erlang 
B formula in (3.2), and the only difference between the limits of Dp and B is the factor 1 + F(l), 
which incorporates all information about the admission policy. 

Let D F (s,\) = D F (s,X + n) and D§(s,X) = D^(s,X + n) with O as in (2.13). Hence, D F (s,X) 
and Dp(s, A) are the stationary probability that an arriving customer finds all servers occupied, 
and the stationary probability that an arriving customers is rejected, respectively, in the system 
with retrials using the RTA approximation. 



Therefore, 



and 



S 



Theorem 4 (QED limits with retrials). Under condition (2.5), for A = s — ^fyfs, with 7 G (—00, 
fixed, and with Vt defined as in (2.13), 

lim = a (3.14) 

s— >oo s 



with a the unique positive solution of the equation 

= f(l ~ a ) 
$(7 - a) ' 



Furthermore, 



(3.15) 



lim ^- s D F (s,s- 1 V~s) = (l + F(l)) ^ a) (3.16) 

lim yfiDgis, s-jy/s) = (3.17) 
s-s>oo $(7 — a) 

The proof of Theorem 4 is presented in Section 7. Theorem 4 shows that the additional load 
due to retrials $7, for a system with many servers, is of the order yfs. In particular, as the number 
of servers grows large, O is well approximated by ay/s, where a is a constant that no longer depends 
on s. This also means that for the overall retrial system the arrival rate A + f2 is approximately 
s-(j-a)y/s. 



4 Dimensioning problems 

First consider the situation without retrials, and the problem of finding the arrival rate A such 
that the probability Dp(s,X) to find all servers occupied or the probability Dp(s,X) that service 
is denied altogether has a prescribed value. Here the number of servers and the admission policy, 
embodied by F, are assumed to be given. 
Problem 1 For fixed s, e, find 7 such that 

y/s~D(s, s - 7^) = e with D = D F or D§. (4.1) 

In Section 6 it will be shown that, under condition (2.5), Dp(s, s — jy/s) decreases strictly from 
1 at 7 = 7p to at 7 = y/s. Unfortunately, such a result does not hold for Dp: There are policies 
F satisfying condition (2.5) such that Dp(s,s — jy/s) is not monotonic as a function of 7, see 
Appendix A. These policies are, however, rather rare. Relevant policies for which Dp(s, s — jy/s) 
is monotonic include 

p k =pe (0,1), k^s, (4.2) 

and, for some N ^ s, 

( 1, s ^ k < N, ,. 

*={o, k>N+i. (4 - 3) 

When Dp(s, s — Jy/s) is monotonic, it decreases from 1 — P at 7 = jp to at 7 = yfs. 
We have the following result: 
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Theorem 5 (Unique solutions). Under condition (2.5), 

(i) Equation (4.1) with D = Dp has a unique solution 7 = 7 Si _f(£) for any e S (0, yfs). 
Assuming further that Dp(s, s — jy/s) is monotonic in 7, 

(ii) Equation (4.1) with D = Dp has a unique solution 7 = 7^ F (e) for any e £ (0, (1 — P)\fs). 

We next consider Problem 1 in the QED regime, and first introduce some definitions. Let for 

7 € (jp,y/s) 

9s,f(7) ■= VsD F (s, s - jx/s), g*p{i) := \/sD§(s, s - 7x /s). (4.4) 
Furthermore, define for 7 € R, 

.00(7) = (4.5) 

and define 700, f(e) and 7^(e) as the solutions of 

(l + F(l)) 5oo ( 7 ) =e and 500(7) = e, (4.6) 

respectively. It is well-known that 500(7) strictly decreases from +00 at 7 = —00 to at 7 = 00, 
and so both equations in (4.6) have unique solutions when e > 0. 

In Theorem 7 below, we give a limit result for 7 Sj _p(e) and r y?p(s) as s — > 00 that involves 
7oo,f(e) and 7^ F (e), respectively. For this result, the following observations are made. From (2.8) 
and (2.10) we get 

^(7) = V~ sB (s, s - lV~s) ^ R( l + F{l ~^lf ) 1 ry (4.7) 

l + B(s,S-^y/s)F(l-j/y/s) 

and 

l + B(s,s-7Vs)-F(l-7/Vs) 

Now we have from [9, Theorem 14], 

v^B(s, s - 7V^) = 5oo(7) + ^^00(7) + 0{s~ l ), (4.9) 

V s 

where 

hoc (7) = " I (7 3 + (7 2 + 2) 5 oo (7)) 9oo (7) , (4.10) 
and the O in (4.9) holds uniformly in any bounded set of 7's. Using (4.9), together with 

F(l -j.) = F(l) - -^F'(l) + 0(s^), (4.11) 

which holds because of assumption (2.5), the following result is established upon computation. 
Lemma 6. Under condition (2.5), 

9s,f(7) = (1 + ^(l))5oo(7) + + O(s^), (4.12) 

9?,fW = 2oo(7) + 4^(7) + Ois- 1 ), (4.13) 
V s 
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where 



W(7) = (1 + ni))M7) - (7^(1) + (1 + F(l))F(l)) goo (l), (4-14) 
^,f(7) = M7) " (7 + 5oo(7))3oo(7)^(l)- (4.15) 

The O in (4.12) and (4.13) holds uniformly in any bounded set of'f's. 

Theorem 7 (Asymptotic dimensioning without retrials). Under condition (2.5), 

1s,f{ £ ) = 7oo,F(e) + A=Voo,f( £ ) + Oi 3 " 1 ), (4-16) 
Vs 



with 



7?oc,f(e = -77— TvTvTTT TTTT' (4-18) 

"- F<E) = - ^(TgM) ■ (4 ' 19> 

Theorem 7 gives the limits of 7 S] ^(e) and 7^ F (g) as s — >• 00, together with first-order correc- 
tions 77oo,iK e ) an d f( £ ) * na * are simple functions of the limits 7oo,f(£) and 7^(e). It is here 
instrumental to note that 

<4(7) = -500(7) (7 + good))- (4.20) 

The proof of Theorem 7 will be given in Section 8. 

We next discuss the numerical experiments that we conducted to illustrate the analytical results. 
Remember that the objective of Problem 1 is to determine maximal sustainable load A such that 
the rejection probability Dp(s,X) is below a threshold e s := e/yfs. In Section 1 we have denoted 
the true maximal load by A op t, and we have explained the concepts of square-root staffing and 
asymptotic dimensioning, in order to obtain accurate estimates of A op t that are asymptotically 
sharp in the QED regime. 

The conventional square-root staffing rule is to use the QED approximation ^Dp(s, A) ~ 
.D*(7), obtain the solution to D*(-y) = e s , say 7*, and then prescribe the load as A* = s— 7*\A- The- 
orem 7 allows for refined-square root staffing based on a better QED approximation ^Dp(s, A) ~ 
D,(^) and the solution 7. to .0.(7) = e s . 

For the asymptotic dimensioning sketched above and in Section 1, applied to Problem 1, we 
identify the following key functions and parameters: 

0.(7) = (4-21) 
1 

7) + " 

7* = 7^), (4.23) 
7. = 7* + -^<f(£), (4-24) 



£.( 7 ) = Z)*( 7 ) + -W£ ( 7 ), (4-22 
V s 
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and using the square-root rule A = s — "fy/s, 



\± = s-^y/s, (4-25) 

A. = s — 7.\/s = A* + r., (4.26) 

with 

r. = 44^. (4.27) 

Table 1 presents results for the admission policy Pk = P = 0.1, k ^ s, and for Problem 1 with 
s = 100 servers. Note that |A op t — A,| is always less than 0.1, and how the refinements r, lead to 
much sharper estimates of the true optimal values. 



£ 


A pt 


A* 


A. 


r. 






0.010 


75. 


,324 


72.836 


75.409 


2.573 


0.004 


0.010 


0.020 


77 


.554 


75.504 


77.621 


2.117 


0.011 


0.020 


0.030 


78 


.996 


77.201 


79.053 


1.852 


0.018 


0.030 


0.040 


80 


.096 


78.479 


80.146 


1.667 


0.026 


0.041 


0.050 


80 


,999 


79.519 


81.045 


1.525 


0.034 


0.051 


0.060 


81 


.774 


80.405 


81.816 


1.411 


0.043 


0.061 


0.070 


82 


,458 


81.181 


82.497 


1.315 


0.052 


0.071 


0.080 


83 


,073 


81.876 


83.110 


1.234 


0.061 


0.081 


0.090 


83 


,636 


82.507 


83.671 


1.164 


0.071 


0.091 


0.100 


84 


,157 


83.088 


84.190 


1.102 


0.080 


0.101 



Table 1: Results for the admission policy Pk = P = 0.1, k ^ s, and for Problem 1 with D = Dp 
and s = 100 servers. 

Table 2 presents results for the same situation, except with Pk = P = 0.5, k ^ s. This control 
policy is thus more lenient, in the sense that when all servers are occupied, 50% (instead of 10%) 
of the arrivals is admitted. Note that despite the fact that the refinements r, are larger, and 
I A op t — A, | is always less than 1, the estimates are less accurate. Based on this example, and many 
other numerical experiments not reported here, we conclude that square-root staffing becomes less 
accurate for systems with admission control policies that allow more customers to enter. Such 
systems really benefit from the refined staffing rules. 

Let us now turn to the same dimensioning problem, but then for the system with retrials. 
Problem 2 For fixed s,e, find A such that 

y/sD§ (s, A + n) = e (4.28) 

with f2 defined as in (2.13). 

We only consider the case D§ in (4.28), but results for D F can be obtained in a similar manner. 
Let us first bring the generalized Cohen equation (2.13) in a form that is amenable for analysis in 
the QED regime. Write $7 = a-^/s so that (2.13) becomes 

ay/a = (s - (7 - a)y/a) D§(s, s - (7 - a)y/s~). (4.29) 
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g A opt X* A. r, t/sD§(s, Ah.) y/s£>ff(s,A.) 



0.010 75.910 72.836 76.298 3.462 0.003 0.011 

0.020 78.162 75.504 78.510 3.006 0.009 0.022 

0.030 79.619 77.201 79.942 2.741 0.015 0.033 

0.040 80.730 78.479 81.035 2.556 0.022 0.043 

0.050 81.642 79.519 81.933 2.414 0.029 0.054 

0.060 82.425 80.405 82.705 2.300 0.037 0.064 

0.070 83.116 81.181 83.386 2.204 0.045 0.074 

0.080 83.738 81.876 83.999 2.123 0.053 0.084 

0.090 84.307 82.507 84.560 2.053 0.061 0.095 

0.100 84.832 83.088 85.078 1.991 0.070 0.105 



Table 2: Results for the admission policy Pk = P = 0.5, k ^ s, and for Problem 1 with 
servers. 

With ff F defined as 



we can write (4.29) concisely as 

a = SsAl ~ a) 

in which 7 is given and a is to be solved. 

We thus have D F (s, A + Q) = e/yfs if and only if 

9?,f{i ~ a*,F(n)) = £, 

where a = a Sj _p(7) solves for 7 > the equation (4.31). 

Using f^ F (S) = (1 — 6/y/s) g^ F (S), the equation in (4.32) takes the form 



1 - (7 - a)/JS 

where a = a s f{i)- This gives 

Q= 1 -^ £ 

and therefore 

7-e 

7 — a 



1-e/y/s 

Hence, we can solve Problem 2 by first solving 5 = S^ F (e) from 
see (4.32), and then setting (7 — e)/(l — £y/s) = 5, i.e., 
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As to (4.36) we operate under the assumption of monotonicity of Dp as made in the beginning of 
this section. 

Denote by 5 = <5oo(e) the solution of 

9oo(5)=e. (4.38) 
Observe that <5oo(e) = 7m( £ )> see (4-6). Then we have by Theorem 7, (4.34) and (4.35) that 

^=^£-^ ,' , x t-Ts +0{s . 4.39 
Using this in (4.37), we arrive at the following result. 

Theorem 8 (Asymptotic dimensioning with retrials). Under condition (2.5), and assuming that 
Dp(s, s — 7\/s) is monotonia in 7, 

^ F {e) =6 + 600(6) + -U* F (e) + 0(s~ 1 ), ( 4 - 4 0) 

0£A £ ) = - S oo(s)e ' ■ 4.41 

For the asymptotic dimensioning scheme, applied to Problem 1 with retrials, we identify the 
following key functions and parameters: 

7* = 8 + 5^(8), (4.42) 

7.=7* + 4=0£,f( £ ). ( 4 - 43 ) 



and using the square-root rule A = s — "fy/s, 



\± = s - 7*a/s ; (4.44) 
A. = s-7.Vs = A* + r„ (4.45) 



with 



, > , ^,F(7*-g) ... 

r. = (7* - e)e + —+- -. (4.46) 

9'oo(l* -e) 

Table 3 presents results for the admission policy Pk = P = 0.1, k ^ s, and for Problem 1 with 
retrials and s = 100 servers. Table 4 displays the results for p = 0.5. 



5 Bistability 

In Section 4 we have considered dimensioning problems that were formulated directly in terms of 
the probabilities Dp and Dp, and we have seen that these problems admit generically one solution 
when the prescribed value for Dp or Dp is in the appropriate range. The situation is different when 
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£ 


Aopt 


A* 


A. 


r. 


y/sD§(s,K) 


v/i£$(s,A.) 


0. 


010 


75.249 


72.736 


75.336 


2 


.600 


0.004 


0.010 


0. 


020 


77.399 


75.304 


77.470 


2 


.166 


0.010 


0.020 


0. 


030 


78.759 


76.901 


78.822 


1 


.921 


0.018 


0.031 


0. 


040 


79.775 


78.079 


79.832 


i 
i 


7^ 
. t oo 


0.025 


0.041 


0, 


050 


80.594 


79.019 


80.647 


1 


.628 


0.034 


0.051 


0. 


060 


81.283 


79.805 


81.333 


1 


.528 


0.042 


0.061 


0. 


070 


81.880 


80.481 


81.929 


1 


.447 


0.051 


0.071 


0. 


080 


82.409 


81.076 


82.455 


1 


.379 


0.059 


0.081 


0. 


090 


82.884 


81.607 


82.929 


1 


.321 


0.068 


0.091 


0. 


100 


83.315 


82.088 


83.359 


1 


.271 


0.077 


0.101 



Table 3: Results for the admission policy Pk = p = 0.1, k ^ s, and for Problem 2 with retrials and 
s = 100 servers. 





£ 


A pt 


A* 


A. 


r. 









.010 


75.834 


72.736 


76.225 


3.489 


0.003 


0.011 





.020 


78.006 


75.304 


78.359 


3.055 


0.009 


0.022 





.030 


79.380 


76.901 


79.711 


2.810 


0.015 


0.033 





.040 


80.407 


78.079 


80.721 


2.642 


0.021 


0.043 





.050 


81.234 


79.019 


81.536 


2.516 


0.028 


0.054 





.060 


81.930 


79.805 


82.222 


2.417 


0.036 


0.064 





.070 


82.534 


80.481 


82.817 


2.336 


0.043 


0.074 





.080 


83.068 


81.076 


83.344 


2.268 


0.051 


0.085 





.090 


83.548 


81.607 


83.817 


2.210 


0.059 


0.095 





,100 


83.984 


82.088 


84.248 


2.160 


0.067 


0.105 



Table 4: Results for the admission policy Pk = P = 0.5, k s, and for Problem 2 with retrials and 
s = 100 servers. 

dimensioning problems of a more complicated nature are considered. In this section we consider 
carried traffic quantities A(l — D(s,X)) with D = Dp or Dp, without and with retrials, and we 
will see that the corresponding dimensioning problems may have two solutions. One could refer to 
this situation as bistability, which has also been observed in certain loss networks with alternative 
two-link routing [7]. The bistability stresses the fact that sample paths of the Markov process tend 
to be concentrated around two relatively stable points of the system. 

Note that A(l — Dp(s,X)) is the rate of arrivals that enter the system in one attempt, and 
A(l — Dp(s,X)) is the rate of arrivals that pass the system, possibly after having to wait in the 
queue that arises when all servers are occupied. 

We first consider the system without retrials. 
Problem 3 For fixed s, e, find A such that 

A(l - D(s, A)) = £ with D = Dp or D§. (5.1) 

In Section 6 we prove the following result. 
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Theorem 9 ((Non-)uniqueness of solutions). Under condition (2.5), 

(i) Equation (5.1) with D = Dp has at least two solutions A E (0,Ap) when e > is sufficiently 
small. 

(ii) Equation (5.1) with D = Dp has a unique solution A G (0, Ap) u>/ien e G (0,s). 

Theorem 9(i) is intuitively clear since A(l — Dp(s,\)) is positive for all A 6 (0, Ap) while it 
vanishes when A | or A | Ap. The result of Theorem 9 (ii) is less obvious and depends on 
monotonicity of A(l — Dp(s, A)) as a function of A £ (0, Ap). For the case pk = p £ (0, 1), k ^ s, 
it is shown in Appendix A that A(l — Dp(s, A)) is strictly concave, and so (5.1) with D = Dp has 
exactly two solutions when e > is sufficiently small. 

We next consider carried traffic quantities in the case of retrials. For convenience, we shall 
restrict ourselves to the choice D = Dp in (5.1). 
Problem 4 For fixed s,e, find A such that 

A(l A + «))=£ (5.2) 

with f2 defined as in (2.13). 

Observe that A(l — Dp(s, A + O)) = e if and only if 

(s - 7V ^) (l - -j= g* F {i - a s ,p( 7 ))) = e. (5.3) 

Using f^ F {5) = (1 — 5/yfs) gfp{b~) and (4.31) in (5.3), we can write equation (5.2) as 

(.-■^(l-l ° ) = ^ s ^jL =e , (5.4) 

V ^/s 1 - (7 - a)/^sJ V s ~7 + a 

where a = a Sj p(7). It is now not possible to simply eliminate a, and a study of the function 

^(7) = ^7 7)2 ' (5.5) 

is required. 

We now give a detailed presentation of what can be achieved analytically for the case that F = 
(Pk = for all k ^ s). It is a challenging problem to explore in what respect this detailed analytic 
result can be extended to more general admission policies. We thus consider the carried-traffic 
problem in (5.3) for the case that F = and study the function 

L s ( 7 )=L s ,p(7)= r (y/ ~ S ~J )2 , r , 0< 7 <v^ (5.6) 
Vs-7 + a s ( 7 ) 

where a s (7) is the solution of the Cohen equation a = f s (^y — a) with f s = f s ,F=o- Both f s and a s 
have been studied in great detail in [2] . Using the results of [2] and extensions thereof, the following 
is shown in Appendix B. 

Theorem 10. There holds 

L s (7)e(0,v^- 7 ), 7 6(0,VS), (5.7) 
L.(i) = 7 5(l + 0( 7 v/i)), 710, (5.8) 

L s ( 7 ) = (v^-7)(l + 0(^(l--^) S )), 7tv^- (5.9) 
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Furthermore, L s (j) is unimodal on (0,y/s), and the maximum of L s {^f) is assumed at the unique 
solution 7 = 7s of the equation 

7«.(7) = i(l-^), (5-10) 



and 



Finally, 



L s {ls 



1, 



— (Vs-%)- 



(5.11) 



1 \l/2 



7»>/.0)=%*>G + Ti) 



1 



s > 2, 



.2 ' 16s/ 

and 7 S increases in s from 1/2 ai s = 1 to 1.034113461 at s = oo, with % ~ 1 lo/iera s = 550 
From the detailed information provided by Theorem 10, it is seen that the equation 

Cv^-7) 2 e 



(5.12) 



^/i - 7 + a ^/i' 



(5.13) 



see (5.4), has two, one or zero solutions according as e <, = or > ■ v /iLs(7s) with L s (^ s ) given in 
(5.11). Furthermore, it is seen that ^/sL s (^ s ) = s + 0(yfs), and that the highest carried-traffic 
numbers occur in a 7-region somewhat away from, but relatively close to, 7 = 0. 

Theorem 11 (Bistability). Problem 4 has two, one or zero solutions according as e <, = or > 
y/s L s (7s) with L s (j s ) given in (5.11). 

In Figure 1 we display 4= L s (5^), < 5 < 1, for s = 1, 5, 10, 50, 100. It is observed that the 
graphs approximate the graph of the function 1 — 5, < 5 < 1, when s gets large. 



1-0 r 




Figure 1: Plot of -4= L s (5y/s), < 5 < 1, for the cases s = 1, 5, 10, 50, 100 , with maximizer 7, 
decreasing with s and maximum value -7= L s (7s) increasing with s. 



17 



6 Proof of results on dimensioning 

In this section we drop condition (2.5) so that both P = 1 and F(Xp — 0) < oo are allowed. 
The equations in (4.1) and in (5.1) take the form 

ail) = e, 9 = 9s,F or gf F , (6.1) 

and 

7 + /( 7 ) = ^£ f = f OTf R (6.2) 

respectively, where gf F and as in (4.4) and (4.30), and we define for 7 £ ( 7 p, y^i], 

/ s ,f(7) == (l " ^(7) = V^(l - £> F ( S , - ls Ts). (6.3) 

To get insight and pertinent results about when the equations in (6.1) and (6.2) have (unique) 
solutions and how to find these solutions, we relate the functions f s ,p and f^ F to the functions f s 
and g s given by 

/ s ( 7 ) = (l-^)^(7) = V^(l-^)s(a,s-7V^), 7 < v 7 ^- (6.4) 

The latter functions have been studied in considerable detail in [2], in particular with respect to 
monotonicity properties. Some of these properties are collected in the beginning of Appendix A. 

Lemma 12. For 7 £ {jp,y/s\, 

1 +(l-^)^(7) 

i + _/ s ( 7 )i?- s ( 7 ) 

and 

1-4=^.(7) 

/£f(7) = / s (7) , (6.6) 

1+ / s ( 7 )iT s ( 7 ) 
Vs 

with 



H s (j) = ^p s ----p s+n (l--t=) . (6.7) 

n=0 

Proof. This follows in a straightforward manner from (2.8) and (2.10) and the definition of H s , 
where we note that F(l - ^) = (1 - # s ( 7 ). □ 

Lemma 12 shows that / Si p and f^ F factorize into an admission policy independent part / s ( 7 ) 
and a part comprising the admission policy via H s . By dividing either side of (6.5) and (6.6) by 
(1 — "f/y/s), it is seen that a similar result as Lemma 12 holds for g s ^ F and g^ F , with the same 
factors comprising the admission policy as in (6.5) and (6.6). 

The following result gives a global picture for f s> p and f^ F in terms of inequalities. We observe 
that F(X P - 0) = 00 H s (^p + 0) = 00. 
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Proposition 13. (i) For 7p < 7 < yfs, 

max{0, -7} < /*,( 7 ) < f s {l) ^ fs, F (l) < ™in{^ - 7, ^ffZ) }- (6.8) 

(ii) There is equality in the first inequality in (6.8) if and only if 7 £ (0, -^/i) and pj. = 1 /or a// 
k ^ s, and in that case 

fs 

H{-/) = —, 7 p = 0<7<Vs. (6.9) 

7 

(mj There is equality in the second inequality in (6.8) for any 7 £ (7P, -^/i) i/ and on/y if Pk = 
for all k ^ s. There is equality in the third inequality in (6.8) /or any 7 £ (7p,y / s) i/ and 
only if pk = /or all k ^ s. 

(iv) There is equality in the fourth inequality in (6.8) if and only if '7 £ (0, -^/s) and p^ = 1 for all 
k ^ s. 

(uj For7 = ^/s, 

/?f(7) = fsiri) = fsAi) = °- ( 6 - 10 ) 

( vi ) /^f(7p + 0) = -7p i/ and onZy ?/ i? s (7p + 0) =00. 

(vii) f s ,F(jp + 0) = \fs — 7p i/ and onfa/ if H s (-f P + 0) = 00. 

The proof of Proposition 13 is given in Appendix A and uses the representations (6.5) and (6.6). 



v 3.0 




\>V 2.5 




2.0 




^^O^W \v 1-5 





-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 



Figure 2: Picture illustrating / s (? , f Sk p and / s with bounds as in Proposition 13 for the case that 
< P < 1 and H s {jp + 0) = oo (actual choice: s = 1, Pk = p = 1/3, k ^ s) 

In Figure 2 we show /fp, f s and / s .p, together with the bounds in (6.8) for the case that 
— oo < 7p < and H(^p + 0) = oo (actual choice: s = 1 and pk = p = 1/3, k ^ s). 

In general, monotonicity properties for the functions f s ,F, gfp and of the functions 7 + / s ,p(7) 
and 7 + //p(7), see (6.1) and (6.2) are not easy to establish or manifestly not true. We now present 
some results, positive and negative, on monotonicity and convexity of the functions in (6.3) and 
(4.30). 
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Proposition 14. (i) f s p and g S) p are strictly decreasing. 

(ii) 7 + //V(t)j 7p < 7 ^ v/s, is strictly increasing. 

The proof of Proposition 14 is given in Appendix A. An interesting example is obtained when 
we choose 

p fc =p 6(0,1), fc^s. (6.11) 
Then P = p, jp = — yfs £ (— oo, 0), and 



p 

oo 



.12) 



Furthermore, H s {^p + 0) = oo, and we compute from (6.5), (6.6) and (6.12) 

Js,f(V v—, TTT' Js,fK1) ~ p— — ■ (o-l^J 

There is for this case the following result. 

Proposition 15. Assume ft=pG (0,1), k ^ s, so that f s p and f^ F are given by (6.13), and let 
g S: p and g R F be the functions associated to f s ^p and f^ F according to (6.3) and (4.30) , respectively. 
Then, 

(i) fs,F an d f^p are strictly decreasing and convex in 7p < 7 ^ yfs, with 

IsAlP + 0) = - 7 P = — , f? F (-yp + 0) = -7P = — A (6-14) 

p p 

and fsAV^) = fSpiVs) = 0. 
(H) 9s, F md g^ F are strictly decreasing in "fp < 7 ^ ^/s, with 

9s,f{ip + 0) = v 7 ^, ^p(7p + 0) = (1 - p) y^, (6.15) 

and g s ,F(Vs) = gf iF (\/s) = °- 
fm,) 7 + f^pij) is strictly increasing and convex in 7p < 7 ^ u>i£/i 

7 + / s R p(7) , =0, 7 + / s fl p(7) _ =v^- (6.16) 



7=7p+0 



7=Vs 



(wj 7 + fs,F(l) is non-monotonic and strictly convex in jp < 7 ^ ^/i, u;zf/j 

7 + /-,f(7) =7 + /,,p(7) = v^- (6.17) 



7=7P+0 



The proof of Proposition 15 is given in Appendix A. 
We finally have the following result. 

Proposition 16. Assume that pk is chosen as in (4.3). Then gf F {^) is strictly decreasing in 

7 e (7p, V~s). 

The proof of Proposition 16 is given in the Appendix A. 

We now show how the various results given in this section can be used to prove Theorem 5 and 
Theorem 9. 
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6.1 Proof of Theorem 5 

Under condition (2.5) we have F(Xp — 0) = H s {^p + 0) = oo. It follows from (6.3) and (6.5) that 

1-7P/V S 1 - 1P/V S fs{lp)/V s 
In a similar fashion it follows from (4.4), (4.30) and (6.6), using 7p = -(1 - P)P" 1 1 /«, that 

g* F (lP + 0) = (l-P)V~s. (6.19) 

Furthermore, 

D F (s, A) = v^Ht), D§(s, A) = xTsgf^) (6-20) 

when A = s — jy/s while g s ^F{\/s) = g^p{y/s) = as B(s, 0) = 0. Then Theorem 5(i) follows from 
Proposition 14(i) while Theorem 5(h) holds by the monotonicity assumption made in the discussion 
preceding Theorem 5. 

6.2 Proof of Theorem 9 

The function A(l — Dp(s,X)) depends continuously on A G (0, Ap), is positive for A G (0, Ap), and 
satisfies 

lim A(l - D F (s, A)) = = lim A(l - D F (s, A)), (6.21) 

A4-0 A'fA p 

see (6.18) and (6.20). This yields assertion (i). 

A computation, using (4.30) and (6.3), shows that 

A(l- J D( S ,A)) = S - v / i(7 + /(7)), (6-22) 

when A = s — J\fs. Hence, by Proposition 14(h), we have that A(l — Dp(s, A)) is strictly increasing 
in A G (0, Ap). Furthermore, 

limA(l - £>p(s,A)) = 0, lim A(l - Dp(s, A)) = s (6.23) 

A4-0 A'fAp 

by (6.18) and (6.20). This yields assertion (ii). 

Let us conclude this section with the following two observations. Concavity of A(l — D(s,X)) 
when pk = p G (0, 1), k ^ s, follows from (6.22) and Proposition 15(h). Monotonicity of D§(s, A) 
when pk is given by (6.3) follows from (6.20) and Proposition 16. 

7 Proof of Theorem 4 

We prove Theorem 4 under condition (2.5), where we assume P G (0, 1) (the case P = requires 
only minor modification of the analysis below). By Theorem 1, the generalized Cohen equation 
(2.13), written in QED coordinates as in (4.31), has a unique solution a s f(t) f° r an Y 7 £ (0, yfs)- 

Theorem 17. For any 7 > 0, 

^00(7) > o,s,f{i) = 000(7) + °(~f=)' s ^ 00, (7.1) 
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where the O holds uniformly in any compact set of 7 G (0, 00), and where 000(7) is the unique 
solution of 

• = w-) = |^. < 7 - 2 > 

see [2, Section 3.3] . 

Proof We have f^ F (5) ^ / s (o") < foo($) and = /oo(^) + 0(775) uniformly in any compact 

set of 5 G E, see (6.6), Proposition 13(i) and [2], Theorem 13 and 15. Let 7 > 0. We prove below 
that 0^(7) 000(7)- Therefore, there is an M > such that 

/oo( 7 - a) - M= < 75,(7 - a) < /oo(7 - a) (7.3) 

holds for a G [0, 000(7)] and all s ^ 1. Since a + /oo(7 — a) is convex in a and f'^S) > —1, 5 £ R, 
see [2, (3.30)], it follows from the mean value theorem that there is an n > such that 

Zoo (7 - a ) - a ^ ??( fl oo(7) - a) 5 a G [0,aoo(7)] > (7-4) 
where it is observed that /oo(7 — a) — a = at a = 000(7)- Then 

= &(l ~ a s,F(l)) ~ o s ,f(7) ^ /oo(7- a s ,jr(7)) - «s,f(7) 7= 

M 

^ ??(Ooo(7) - o,s,F(l)) ~ —Fi ( 7 - 5 ) 

Vs 

and so 

M 

a s ,F(7) > 000(7) /=■ (7-6) 

In the above argument, M and 7/ can be chosen independently of 7 in any compact subset of (0, 00), 
and so the result follows. 

We still have to show that 0^(7) ^ 000(7)- We have from f^ F (5) ^ fs(8) that 

/oo(7 " a a ,jr( 7 )) > /^(7 - a S!J p( 7 )) = 0^(7). (7.7) 

Therefore, 

7 - «oo(7) + /oo(7 - ^00(7)) =7^7- a s,F(l) + food ~ a s,F(l))- (7.8) 

Now 5 + foo(o~) is increasing, see [2, Eq. (4.8)], and so it follows from (7.8) that 7 — 000(7) ^ 
7 - a SjF (l), i-e-, 0^(7) < 000(7)- □ 

The limiting behavior of Dp(s, A + O) and D F (s, A + Q) is now easily found. We have, see (6.3) 
and (4.12), for 7 > 

—= D F (s, s - (7 - a SiF 7)) = — - — -—— 

V s fs,F{J ~ o s ,f(7J) 

(l + ^(0) ^lo^)) )" 1 = (( 1 + g (°))°°o(7))" 1 , --^oo. (7.9) 
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Similarly, 



1 - (7 ~ Qs,f(7))/\/« 
fR 



-= D F R (s, s - (7 - 0^(7))) 

V s /"f(7-v(7)) 
= -^ rT (l-(7-a S)J ,(7))/%/i)^a- 1 ( 7 ), 00, (7.10) 

where it also has been used that a = a s> F(j) satisfies (4.31). 

8 Proof of Theorem 7 

We give the proof of (4.16), (4.18) in detail; the proof of (4.17) and (4.19) being quite similar. Take 
any 71 , 72 £ K with 

7i < 7oo,F(e) < 72- (8.1) 
From strict decreasingness of 500(7) and the definition of 7o0) j<'(e), we have 

(l + F(l)) 5oo ( 72 ) < (l + F(l)) 5oo ( 7oo , F (e)) = e< (l + F(l)) 5oo ( 7 i). (8.2) 

Because 5^,^(7) — >• (1 + ^(l))5oo(7) as s — )• 00, we have from (8.2) that 

9s,f{12) < e < g s ,F(ji) (8.3) 

when s is large. By monotonicity of g s ,F, this implies that j s> f(s) G [71,72] for large s. Since 71 
and 72 in (8.1) are arbitrary, it follows that 7 Sj i?(e) — > 7oo,f(^) as s — > 00. In particular, j s ,f{£) is 
bounded in s ^ 1. 

We next show from the weakened form 

^(7) = (1 + ^(l))5oo(7) + 0(s^ 2 ) (8.4) 

of (4.12) that 

1,,f(e) = 7oo,f(£)+0(s- 1/2 ). (8.5) 
To that end, we consider (8.4) with 7 = 7 S) ir(e) and write 

5oo(7s,f(£)) = 5oo(7oo,F(e)) + {lsA £ ) ~ looA e ))9AlooA e )) + 0{{-is,F{e) - 7oo,F(e)) 2 )- (8.6) 
Using this, together with 

g s ,F{lsA £ )) = £ = (! + i ? (l))5oo(7oo,F(e)), (8.7) 

in (8.4) we get 

e = £ + (1 + F(l))( % , F (e) ~ ^A^WooilooA^)) + 0((%A £ ) ~ 7oo,F(e)) 2 ) + 0(s^ 2 ). (8.8) 

It is well-known that g'Al) ls negative and bounded away from when 7 is in a bounded set. From 
(8.8) we therefore get (8.5). 
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We finally show (4.16), and for this we repeat the argument for showing (8.5), but now using 
the full strength of (4.12) with 7 = 7 Sj i?(e). Using (8.5) in (8.6) yields 

9oohs,F(e)) = 5oo(7oo,F(e)) + (rfa,F( £ ) - 7oo,F(e))5io(7oo,F(e)) + O(s^). (8.9) 
Furthermore, again by (8.5), 

W(7s,He)) = /ioo,f(7oo,f(£)) + 0(s~ 1/2 ). (8.10) 
When we use this in (4.12) with 7 = 7 S) f(e) together with (8.7), we obtain 

£ = e + (1 + F(l))( 7 , iF (e) - 7oo,F(e))^(7oo,F(e)) + 0(s~ l ) + 4=W(7oo,He)) + O^ 1 ). (8.11) 

From this (4.16) and (4.18) follow at once. 
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A Remaining proofs, except proof of Theorem 10 

We start by recalling some basic properties, shown in [2], of the functions f s and g s given in 
(6.4). We have 

a. is (7) is strictly convex and decreases strictly in —00 < 7 ^ y/s from +00 to 0, and 7 + is (7) 
increases strictly in —00 < 7 ^ y/s from to y/s. Furthermore, 7 + / s (7) = O(-), 7 — > —00. 

b. g s (7) decreases strictly in — 00 < 7 ^ y/s from y/s to 0. 

c. / s (7) and g s {^l) increase in s ^ 1 to foo{l) = #00(7) = <p{l)/^{l) uniformly in any compact set 

of 7 e r. 

Proof of Proposition 13. We have for < 7 < y/s that 

H M < E (l - = A (A.l) 

n=0 v 

with equality if and only if p s + n = 1, n = 0, 1, ... . Hence, 1 — ^ H s (^)/y/s ^ for < 7 < y/s, and 
evidently 1 — ^ H s {^)/y/s ^ for 7p < 7 ^ 0. Therefore, from (6.6), //V(7) ^ 0j 7p < 7 < 
with equality for any 7 if and only if p s + n = 1, n = 0, 1, ... . Next, we write (6.6) for 7p < 7 < y/s 
as 

/&(7)=-7+ r /s(7) ■ ( a - 2 ) 

1+ / s ( 7 )tf s ( 7 ) 

and then it follows from 7 + / s ( 7 ) > that ffp{~f) > —7. This proves the first inequality in (6.8). 
Next, we write (6.6) for 7p < 7 < y/s as 

/&(7) = / s (7) - /.(-T) ^(7) , (A.3) 

V " l + ^/ s ( 7 )# s ( 7 ) 

and it follows from 7 + / s ( 7 ) > that /^p( 7 ) ^ /s(7)> with equality if and only if H s (-y) = if and 
only if p s + n = 0, n = 0, 1, ... . This proves the second inequality in (6.8). 
Next, we write (6.5) for 7p < 7 < y/s as 

1 + f s {j)H s (j)/y/s 
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and it follows from 7 + fs(l) < y/s that f s ,F(l) fs("f), with equality if and only if ^(7) = if 
and only if p s + n = 0, n = 0, 1, ... . This proves the third inequality in (6.8). 
Next, we write (6.5) for 7p < 7 < sfs as 

fsAi) = V? - 7 - f "V/^?v ( A - 5 ) 

and it follows from 7 + / s (7) < \/s that / s ,f(7) < — 7- Since -y/i — 7 > / s (7) for 7 < y^s, we 
have that the function x ^ 1— > (1 + x(l — j/y/s))/(l + x f s {pf)/y/s) is strictly increasing. From 
(6.5) and (A.l) it then follows (with x = H s { 7 ) < yi/ 7 ) that / S)i r( 7 ) < </if s {i)/{l + fs(l)) for 
< 7 < \/s, with equality if and only if p s+n = 0, n = 0, 1, ... . Furthermore, for 7 ^ 0, we have 

r- Vsfsil) 7(\/^-7-/^(7)) . n f A a\ 

V s - 7 1 f - f x = , ; / s < 0, (A.6) 

7 + M7) 7 + M7) 

with equality if and only if 7 = 0. This proves the fourth inequality in (6.8). 

The cases of equality in the inequalities in (6.8) have been indicated already along with their 
proofs, and this settles Proposition 13, (ii)-(iv). 

Proposition 13(v) follows from l(i) and the fact that f s (y/s) = 0. 

Proposition 13(vi)-(vii) follow from the representations (6.6) and (6.5) and the fact that H s (j) 
increases to H s {^p + 0) as 7 decreases to jp by non-negativity of all pk- This completes the proof 
of Proposition 13. 



Note. We have the following consequences of (A. 3), (A. 4) and 7 + / s (7) = O(-), 7 — > —00. 

a- /5,(7) = /*(7) + O( 7 ),7P<7<0, 

b- fsM'y) = fs(Tf) + V^ + O(i), 7 P < 7 < 0, when F £ 0, 

c 9?, F (i),9,,F(i) = y/s + O(-), 7P < 7 < 0. 

These results are in particular relevant when P = so that jp = — oo. 

Proof of Proposition 14. (i) Since / s ,f(7) = (1 — l/y/s) <7s,f(7) ; it is sufficient to show that 
9s,f(i) is strictly decreasing. We have from (A. 5) that 

1 . . 

9s ,F(7) = = 1 ^ • (A.7) 

Now g s (7) strictly decreases in —00 < 7 ^ ^fs from y/s to 0, and f s (7)^(7) is non-negative and 
decreasing in 7p < 7 y/s. It follows that 

(A.8) 

l + _/ s ( 7 )£r s ( 7 ) 
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is non-negative and strictly increasing in 7p < 7 ^ y/s, and the proof is complete. 

(ii) We have that 7 + / s (7) is positive and strictly increasing in —00 < 7 ^ y/s, and / s (7) H a (j) is 

non- negative and decreasing in tp < 7 ^ y/s. It follows that 

\ + (A.9) 

1+ / s ( 7 )tf s ( 7 ) 

is strictly increasing in 7p < 7 ^ y/s. Then it follows from the representation (A. 2) of f^ F that 
7 + ff'pi^f) is strictly increasing in 7p < 7 y/s. 

Proof of Proposition 15. With p^ = p € (0, 1) for k ^ s, we have f s p and f^ F as in (6.13) so 
that, in particular, f^p = (1 —p) f s ,F- 

(i) We have that 7 + /a (7) is positive and strictly increasing in —00 < 7 ^ y/s and that / s (7) 
is positive and strictly decreasing in —00 < 7 ^ y/s. Hence, from (6.13), both f s ,F(j) and f^pij) 
are strictly decreasing in 7p < 7 ^ y/s. We next show (strict) convexity of f a p. We write (6.13) 
for 7p < 7 < y/s as 

f?Ai) = rrrrr-( ^ ; b = n P ^ r = — > °< (a.io) 

l + 6(7 + /«(7)) (I-P)vs 7P 

and we compute for 7p < 7 ^ y/s 

(f R y, M = / /s(7)(l + fc7)~^(7) 
Us ' Fj l7J V (l + K7 + /,(7))) 2 

^(7)(l+ft7)(l+b(7+/ 8 (7)))-2fe(/j(7)(l+ft7)-fc/ 8 (7))(l+/i(7)) 
(l + 6( 7 + / s (7))) 3 

From 



(A.ll) 



£'(7), l + 6 7 >0, / s (7>max{0,- 7 }, - 1< / s ( 7 ) < (A.12) 

for 7p < 7 < y/s, it follows that {ff'p)"{'i) > for 7p < 7 < y/s. Hence, //p(7) is strictly convex 
in 7p < 7 < y/s, and so is / s ,p(7). 

The values of f s ,F, f^p at 7 = 7p + and 7 = y/s, as given in and below (6.14), follow from 
Proposition l(v)-(vii) and the fact that H s (~jp + 0) = 00, see (6.12), in this case. 

(ii) We have for 7P < 7 < yfo from (6.3), (4.30), (6.4) and (6.13) that 

9s Ai) = • (A.i3) 

1 ~P+ —j=P{l + fs{l)) 

Now 5 S (7) is strictly decreasing and positive in —00 < 7 < y/s and 7 + f s in) is strictly increasing 
and positive in —00 < 7 < y/s. Hence, 5^(7) is positive and strictly decreasing in 7p < 7 < y/s, 
and so is g^p(j) = (1 — p) g s ,F(j)- The values assumed by g Sj p and gf F at 7 = 7p + and 7 = y/s 
follow from (i). 

(iii) We have from (i) that 7 + / s F (7) is strictly convex in 7p < 7 < y/s, and (/ s R p)'(7) > — 1 
for 7p < 7 ^ y/s by Proposition 14(h). The values assumed by 7 + f^pij) at 7 = 7p + and 
7 = y/s follow from (i). 
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(iv) Strict convexity of 7 + f St F{j) m 7p < 7 < \fs follows from (i). The values assumed by 
7 + fs,F{j at 7 = 7p + and 7 = \fs also follow from (i). From these values it is seen that 
7 + fs,F{l) is non-monotonic in 7p < 7 < y^. 

This completes the proof of Proposition 15. 



Note. It follows from the expression for (/rp)'(7) in (A. 11) that 

(/^) / (7P + 0) = -^y>-l, (A.14) 

and that 

(/ s ,f)'(7p + 0) = ^— v^r < "I, (A.15) 
where the two inequalities follow from < 5 + f s (5) < ^fs and the definition of 7p, see (A. 10). 

Proof of Proposition 16. With pk = 1, s < k ^ TV, p^ = 0, k > N, we have P = and 
iJ s ( 7 ) = 7~ 1 ^/s(l - (1 - l/\fs) N+1 ). For 7 ^ ^/s it then follows from (6.6) and (4.30) that 

R / \ (1-7/^)^^(7) , A1fi s 

3 s ,f(7) = jf ■ (A.16) 

l + 4/ s (7) ^(1-7/v^r 

V n=0 

A computation shows that 



(i- 7 /Vi)^-i 

= -iV/ s (7)+(l--^)v / ^/;(7)-77|/, 2 (7)E( Ar - n )( 1 -7^) n ' ( AJ7 ) 

n=0 

and all terms on the second line of (A. 17) are negative. Hence, g^ F is strictly decreasing, as required. 
Some examples. Take -£^(7) in (6.7) of the form 

n=0 



with < e < 1, a > 0, < P < 1. Then < q n < 1, <7 n is decreasing in n and lim sup qv/^ n+1 ^ = P. 
When a > 1, we have that H s (^p + 0) < 00. Consequently, in this case 

1 - P < D§{s, X P - 0) < Dp(s, A P - 0) < 1. (A.19) 

When 1 < a < 2 we have that H' s (jp + 0) = —00, and in this case (Dp)'(s, Ap — 0) = —00, showing 
that Df (s,A) is increasing for A close to Ap. An example where F(jj — 0) = H s {^fp + 0) = 00 and 
g^ F (jp + 0) = 00 = -(D§)'(s, Ap - 0) is provided by the choice s = 1 and g n = ^ (~* /2 ) (-|) n , 
n = 0, 1, . . . , yielding ^(7) = i(l + 7 )-V2 wit h P = 2, 7P = -1. 
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Proof of Lemma 2. We have from (6.4) and (6.5) that 



1 n-l/ r\ 1 -7/v / s 

-= D F (s, s - = — — — — 



/ s ( 7 ) 1 + (l- 7 /VS)fl.(7) 

i - 4= 7 #.(7) + 4= KsirtW.irt) + 7) 

V£ y_f 

/.(7) 1 + (1-7/^)^,(7) 



;i - p(7)) 1 s/f + p(7) 4= 1- f ; ; r 

7a (7) V s fs (7) 



with ^(7) given in (3.6). Now, by (6.4), 

l-7/\/^ _ 1 

and 

7 7 



£^(5,5-7^), (A.21) 



(/ a (7) +l) = l- J r + J r B- 1 (s,s-j^) = C-\s, s - 7v^). (A.22) 
fs (7) V« V« 

Then (3.4) follows from (A.20)-(A.22). The proof of (3.5) is similar. 
We observe that ^ p{j) ^ 1, which follows from (A.l). 

Proof of Theorem 1. From [5, Section 3] one can extract a proof of Theorem 1. This proof 
depends on basic properties of birth-death processes. The proof that we give here is presented in 
QED coordinates and uses the analytic properties of f^p, such as given in Propositions 13 and 14. 
The proof is given under general conditions, so that also P = 1 and F(-p — 0),H s (^p + 0) < 00 is 
allowed. 

We distinguish the following cases: 
Case a. P = 0. Then 7p = —00, and both 77 s (7) and 7/^(7) are well-defined and analytic in 
—00 < 7 ^ t/s. Set jp t F = 0. 

Case b. P £ (0,1]. ' Then 7P = -(1 - P)yfs/P £ (-00, 0], and both H a (i) and f^ F (j) are 
well-defined and analytic in 7p < 7 ^ y/s. In this case b, we distinguish the subcases 

bl. H s (-yp + 0) < 00. Then by Abel's theorem, 77 s (7) is continuous in 7p ^ 7 ^ yfs, and so is 
b2. H s (j P + 0) = 00. 

In these cases we have from (6.6) and (A. 2) 

Js,fUP + 0) = -7p H 7 • [A. 26) 

1+ / s ( 7P )tf s (7P + 0) 
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Hence f^ F (jp + 0) > — 7p in case bl while f^ F {^p + 0) = — 7p in case b2. Set 7pp = 
/ s R f(7p + 0)+7p- 

In case P = 0, we have by Proposition 13 and 14 that 

/&(7) > "7, (/^)'(7) > -1, (A.24) 

for 7 < -y/s, and it follows as in the proof of [2], Theorem 10 in Section 4.3 that for any 7 G (0, y/s) 
there is a unique solution a of the equation a = //p(7 — a)- This solution, 0^(7), satisfies 
a s,F(7) - ► +00 as 7 4 0. For if 6 := lim inf 7 ^ o s ,f(7) < °°i we would have 6 = liminf^o f{l ~ 
a s,F(l)) = f(~b), contradicting the first item in (A.24). Similarly, it can be shown, compare the 
beginning of the proof of [2], Theorem 8 in Section 4.8, by considering c := lim supyj.^ 0^(7), 
that a S) ir(7) — > as 7 f \fs- Assume now that P G (0, 1], in which we exclude the case that pk = 1, 
k ^ s. Now (A.24) holds for 7p < 7 < \/s. Let 7 G (7pp, We have 

/ Sj V(7 - 0) > 0, (A.25) 

while 

f^il - (7 - 7p) + 0) = jp, F -IP <1 -IP, (A. 26) 

and so 

f?; F {i-t)<s (A.27) 

when (5 is less than but close to 7 — 7p > 0. By continuity, it follows from (A.25) and (A.27) that 
the equation a = //p(7 — a) has a solution a, and this solution is unique by (A.24). To show that 
this solution, 0^(7), satisfies a S) p(7) — > f^ F {^p + 0) as 7 1 7pp, we need the following lemma. 

Lemma 18. Lei c < < d and let h : (c, d) — )■ R 6e smooth and such that h(a) > 0, /i'(a) < 1 
for c < a < d while h(a) — > d when a\ d. Then for any e, < e < —c, there is a unique solution 
a(e) of the equation a = h(a — e), and a(e) — > d as e J. 0. 

Proof. Let e G (0,-c). Now 

/i(0 - e) - = h(-e) > ; h(a - e) - a -> -e < , a|d + e. (A. 28) 

Hence, by continuity, there is an a G (0, d + e) such that a = h{a — e), and this a is unique by the 
assumption that h'ib) < 1 for b G (c, d). To proceed, we first show that h(a) > a when c < a < d. 
Indeed, when h(a) < a for some a G (c,d), we would have 

h{a\) — h{a) = h(a±) — a\ — (h(a) — a) + a — a\ > a\ — a (A. 29) 

when a\ is sufficiently close to d, contradicting h'(b) < 1 for all b G (c, d). Therefore, h(a) ^ a for 
a G (c, d), and h(a) = a for some a G (c, d) cannot occur either. For otherwise, we would have for 
a\ > a that 

h{a\) — h(a) ^ a± — a, (A. 30) 

contradicting h'(b) < 1 for all b G (c,d). 

We now show that a(e) — > d as e 1 0. We have from a(e) < d + e, see (A. 28), that 
limsup e i a(e) ^ d. Now suppose that b := liminf^o < d. Take e„ > 0, e n — >• such 
that a(e n ) — >• 6 < d. Then 

= h(a(e n ) - e n ) - a(e n ) -> fc(b) - 6, (A.31) 
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contradicting h(b) > b for all b S (c, d). This completes the proof of the lemma. 

Taking in the lemma 

h(a)= f* F {lP,F-a), C :=-V^ + lP,F<a<fK F ( lP + 0)=:d, (A.32) 

it follows that a S: F{l) — > ff'pilP + 0) as 7 | 7p,_f- Finally, 0,5^(7) — > as 7 f y/s can be shown by 
using the same argument as in [2], Theorem 8 in Section 4.8 for showing that 0(7) — > as 7 f \fs. 
This completes the proof. 

B Proof of Theorem 10 

In this appendix we present the proof of Theorem 10 on the function L s in (5.6), given in 
terms of the function a s (j) that solves Cohen's equation a = f s {^ — a). The proofs rely heavily on 
(extensions of) the results in [2]. In particular, we use 

2 7 2x , , / In 7 3 , 



7 0,(7) = 1 - - [1 - -)Y + 4(1 - -J ^= + 0( 7 4 ), 7 I 0, (B.l) 

which is a sharpening of [2, Theorem 5]. This sharpening can be obtained by the method to prove 
[2, Theorem 5] where, as an intermediate step, [2, Theorem 10] should be sharpened to 

/ » w = - s -s-^s + ( 2 -;)i + ( 16 -T)^ +0 (i)- '->-»• (R2 > 

using the methods of [2, Section 4.1]. We shall also use and sharpen [2, Proposition 6], 

2 1 

1 - ^7-7 2 < 7 a s (7) < 1 - -p7, 0<7<\/i. (B.3) 
V s V s 

Proof of < L s (j) < y/s — 7, < 7 < y/s. This follows from the definition in (5.6) and a s (7) > 0. 

Proof of £3(7) = 7s(l + 0(7^/s)), 7 | 0. This follows from the definition in (5.6) and (B.3). 

Proof of L s ( 7 ) = (\/i - 7)(1 + 0(s _1 / 3 e a (l - 7/\/s) s )), 7 t V^- This follows from [2], The- 
orem 8, and Stirling's formula. 

We next show the results on unimodality. 
Proposition 19. We have for < 7 < \fs 

4(7) = O 70,(7) = 5 (1 - 7/v 7 ^)- (B.4) 
We have from the definition of L s in (5.6) 

4( 7 ) = & yfs - 7 + 2a s ( 7 ) + (V5 - 7) 4(7) = 0. (B.5) 
From [2], (4.50) we have 

_ a, W(7 + 1/^) Q < 
1-7/VS-7W) 
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Using this in (B.5) with the facts that 7 > and 1 — j/y/s — 70^(7) > 0, we have for < 7 < \/s 

L'M = & (7 a s (j)) 2 + (1 - 7 /v^)(7^ -|)7 0.(7) " \ U " ll^f 7v^ = 0. (B.7) 
The quadratic in 703(7) occurring in the second proposition in (B.7) has the roots 

7^(7) = -i (1 - 7/\A)(7V^ - |) ± \ (1 - 7/^(7^ + §)■ (B.8) 

Since 70^(7) > 0, only the root in (B.8) with the +-sign needs to be considered. The latter root 
equals | (1 — j/y/s), and this completes the proof. 

To show unimodality of L s , we should consider the function 7a s (7)/(l — y/y/s), < 7 < y/s. 
This function assumes the values 1 and at 7 = 0+ and 7 = y/s — 0, and so, by Proposition 19, 
it is sufficient to show that this function is strictly decreasing in < 7 < y/s. The result we show 
below is somewhat stronger and will also be used in the proof of Proposition 21. 

Proposition 20. 7a s ( 7 )/(l — j/y/s) 2 decreases strictly in < 7 < y/s when s > 1. 

We compute 

((1 - 7/v^r 2 7as(7))' = (1 - 7/v^r 2 [4= (1 " 7/v^r 1 7«s(7) + a, (7) - 7 a' s ( 7 )l • (B.9) 
Using (B.6), we thus see that for < 7 < y/s 
((l- 7 /^)- 2 7as(7))' <0 

(1 - 7 o*(7) " 7/v^) + 1 " 7/>A) - 7(7 + Vv^)(l - 7/v 7 ^) < 

O (l- 7 /v^)(l-7 2 )-7as(7)(l + 7/v^)<0. (B.10) 
It is therefore sufficient to show that for < 7 < -y/i 

7 a s ( 7 )>(l-7 2 ) ?7 7 ^ . (B.H) 
1 + 7/VS 



We compute 



1 -^( 1 -^) = 1 -^-(l-^ 72 + 2fl-I)4 + 0( 7 '), 7 40, (B.12) 



and so, by (B.l), we see that (B.ll) holds for small positive 7. Now suppose that 7, < 7 < 
is such that 

At such a 7 we compute, using (B.6) and (B.12) twice, 

7^(7) (7 + 



(7^(7))' = a s (7) 



1 - l/y/s- 7^(7) 



<,,( 7 )-l±^ fflj(7 ) = -4^L. 

1 - 7/VS Vs 1 + " 
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At the same time, we compute 

\ {1 7) l + 7/^ " V~s (l + 7/^) 2 ' ( ] 

Since s > 1, we have for < 7 < yfs 

1 + -yy/s - 7 2 - 7 3 /v^ > (1 + 7/V«)(l " 7 2 ) = 1 + 7/v^ - 7 2 - 7 3 /v^- (B.16) 
Hence, at a 7 £ (0, y/s) where (B.13) holds, we have 



( 7 a,(7)) / >((l-7 2 ) ?", f 7; ) • (B.17) 



1 + 

This is in particular so for 

70 := inf {0 < 7 < y/s\ (B.13) holds}. (B.18) 

This 70 G (0, y/s) since (B.ll) holds for small positive 7 and since we have assumed that there is a 
7 G (0,y/s) such that (B.13) holds. However, validity of (B.13) and (B.17) for 7 = 70 implies that 

7 a s (7)<(l-7 2 ) |~^ (B.19) 



holds for 7's close to but less than 70. However, (B.ll) holds for 7's close to 0, and so there is, by 
continuity, a 71 < 70 such that (B.13) holds. Contradiction, see (B.18). This proves that (B.ll) 
holds for all 7 G (0, v^)> and the proof is complete. 
Note. We have 703=1(7) = (1 - j) 2 , < 7 ^ 1. 

We have from Proposition 20 that there is a unique root 7 = j s G (0, \/s) of 

7 o s ( 7 ) = |(l-7/v^)- (B.20) 

From the behavior of L s {^f) near 7 = and 7 = y/s, it thus follows that L s is unimodal, with 
unique maximum at 7 = The value of £3(7) at 7 = see (5.11), is easily obtained by inserting 
(B.20) into the definition of L s in (5.6). 

We proceed by showing some properties of %. 



Proposition 21. For s > 1, 



Proof. It follows from a s (j) = f s (j — 0^(7)) that a s (f s (0)) = / s (0). Using the inequality in (B.ll) 
with 7 = f s (0), we get 



(0) = (1 - tf(0)) (B.22) 
1 + 7s(0)/V s 



Working this out (noting that the cubic terms fg (0) cancel), we find 

f?f0)>ifl-^ 



#(0)>-(l-^). (B.23) 
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That is, the value assumed by 7 a s (7)/(l — 7/- v /i) at 7 = / s (0) exceeds 1/2. Since 7 a s (7)/(l — "f/y/s) 
is strictly decreasing in < 7 < \/s, the first inequality in (B.21) follows. Finally, from (B.23) we 
also get 

«°> + 47i) 2 " \ ~ m > °- (B.24) 
and this gives the second inequality in (B.21). □ 
Note. We have % =1 = 1/2. 

Theorem 22. 7 S increases strictly in s ^ 1 /rom 1/2 oi s = 1 to 1.034113461... at s = 00. 
Proof. Let s > t ^ 1. From (B.3) we have for u 1 that 

Pu(7) == i 7 = 1 ~ + Q(7 2 ), 7 10, (B.25) 

1 - 1/v« 



and so 

7« s (7) > W7) (B26) 



1 - 7/Vs 1 - tM 

holds for small positive 7. Now suppose we have a 7 G (0, such that 

p;= _7^ = _W2L (B . 27) 

1-7/V* 1-7/v* 
We have p € (0, 1). We compute, using (B.6) and the definition of p u in (B.25) 

o u (7) V_ _a "(7) 7 + 1^ 1 au(7) 



1-j/y/u 1 - 7/\/u- 70^(7) ( 1 -7/\/^) 2 
" Pui7j 7 (l-7/v^)(l-Pn(7)) 7(l-7/v^)' ^ j 



Hence, at a 7 where (B.27) holds, we have 

0,4(7) V -p 7 + l/y/u p 1 



.1-7/^/u/ 7(1 -p) l — 'y/y/u 7^/n l-7/v^' 
Therefore, at such a 7 we have 

a* (7) V / a t(7) 



u = s,t. (B.29) 



p p + i/y/s _ 7 + i/yt ^ p 



7(1 -P) ^1 -l/y/s l-j/y/i/ 7^\/^-7 V* 
p(#-VVi) / 1 + 7 2 



7 (l- 7 /^(l- 7 /Vt) V 1-P 



+ 1 > (B.30) 



since s > t, p £ (0,1) and 7 G (0, Letting 70 = inf{0 < 7 < >/*| (B.27) holds} G (0,y/i), 
we arrive at a contradiction in the same way as in the proof of (B.ll). Hence, (B.26) holds for all 
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7 £ (0, y/i). Then from decreasingness of 7 o s (7)/(l — 7/ as a function of 7, it is then seen that 

7s > % 

The limit value of -% as s — > 00 follows on considering the equation 7000(7) = 1/2 in which 
000(7) is the unique solution of the Cohen equation in (7.2). We still need a lemma. 

Lemma 23. 7000(7) strictly decreases in < 7 < 00 /rom 1 to 0. 

Proof. Using 

/ 1 x -7^00(7) m Q -n 

°- (7) = l-7aoo(7) ' (R31) 
compare (B.6), it follows for 7 > that 

(7000(7))' < 0^7^00(7) > 1"7 2 - (B.32) 

From [2], Theorem 14, it is seen that 7000(7) > 1 — 7 2 holds for small positive 7. At a 7 > where 
7000(7) = 1 — 7 2 > we compute (7000(7))' = > — 27 = (1 — 7 2 )' and so, by the method of the 
proofs of (B.ll) and (B.26), such a 7 does not exist. Hence, 7000(7) is strictly decreasing in 7 > 0. 
Finally, 7000(7) — > 1 as 7 4- follows from [2], Theorem 14, and 7000(7) — > as 7 — > 00 follows 
from 

7^00(7) = 7/00(7 - 000(7)) < 7/oo(7 - «oo(l)) ->■ 0, 7 ->• 00. (B.33) 
This completes the proof of the lemma. 

We conclude from Lemma 23 that 7000(7) = 1/2 has a unique solution 700, which can be 
determined numerically by a two-stage Newton method, using conveniently 

fM = "/oo(o)(o + /oo(o)) and (B.31) (B.34) 

for solving a = /oo(7 — o) for a = 000(7) an( A 7000(7) = 1/2 for 7, respectively. 

The value of ■js for s = 1 follows from the Note before Theorem 22. □ 
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